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INTRODUCTION 

The  helicopter  rotor  wake  is  extremely  complex  because  of  the  presence  of  multi¬ 
ple  blade  root  and  tip  vortices.  The  interaction  of  these  vortices  causes  a  hiqhly 
distorted  wake  structure,  particularly  at  moderate  helicopter  advance  ratios  (^  =  0 
to  0.2)  and  high  rates  of  descent  of  the  helicopter.  Predictions  of  time-dependent 
and  spatial  velocity  variations  in  this  highly  distorted  flow  field  are,  at  best, 
only  qualitative.  Even  a  basic  understanding  of  the  flow  field  environment  is  lack¬ 
ing,  and  only  a  limited  amount  of  experimental  data  is  available.  Advances  in  dis¬ 
ciplines  that  require  knowledge  of  these  rotor  wake  effects  will  be  limited  until  the 
details  of  the  rotor  wake  environment  are  understood.  These  details  can  only  be 
obtained  through  precise  measurement  of  the  rotor  wake  structure. 

The  available  data  on  rotor  wake  structure  include  measurements  made  with  total- 
pressure  probes,  hot-wire  and  hot-film  anemometers,  photographic  techniques,  local- 
surface-pressure  taps  (from  which  local  in-flow  velocities  can  be  inferred),  and 
advanced  techniques  like  laser  velocimetry.  Each  of  the  data  sets  obtained  with 
these  systems  has  been  valuable  to  the  helicopter  research  community  in  developing  a 
new  appreciation  of  the  complexity  of  the  rotor  wake  structure.  Each,  however,  has 
its  limitations  as  well  as  its  benefits.  No  single  method  can  be  considered  the  best 
for  all  measurement  requirements.  For  instance,  the  laser  velocimeter  (LV)  elimi¬ 
nates  the  influence  of  probe  presence  on  the  velocity  being  measured  and  makes  it 
possible  to  obtain  measurements  between  rotor  blades;  however,  it  is  expensive  and 
can  be  difficult  to  put  into  operation. 

The  LV  is  capable  of  measuring  velocity  time  histories  or  velocity  variation 
with  rotor  azimuth.  However,  this  is  a  difficult  undertaking  because  the  velocity 
measurements  obtained  are  not  made  at  regular  intervals  of  time,  but  occur  only  when 
a  scattering  particle  happens  to  pass  through  the  LV  sample  volume.  The  arrival  of 
these  particles  in  time  can  be  approximated  by  a  Poisson  distribution  function 
(ref.  1).  Successful  experiments  have  been  conducted  (refs.  2  to  4)  using  a  "data 
window"  technique  that  gathers  LV  data  only  when  the  rotor  blade  is  at  a  prescribed 
azimuth.  Another  successful  technique  (refs.  5  and  6)  involves  associating  a  rotor 
blade  azimuth  position  measurement  with  each  sample  velocity  measurement  rather  than 
restricting  measurements  to  a  data  window.  The  latter  technique,  which  is  used  in 
the  Langley  4-  by  7-Meter  Tunnel  to  acquire  LV  data,  is  considered  more  time- 
efficient.  Ibis  report  presents  preliminary  results  of  an  experiment  designed  to  use 
LV  measurements  to  determine  the  effects  of  rotor-fuselage  spacing  and  fuselage  width 
on  rotor  wake  velocity.  Detailed  wake  measurements  were  obtained  to  evaluate  the 
effects  of  these  parameters  on  variations  in  local  time-dependent  velocity  near  the 
rotor  blades. 
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The  axis  systeitr  £he  velocity  measurements  and  measurement  locations 

presented  in  this  r&port- is -presented  in  figure  1(a).  Ihe  values  for  physical  quan¬ 
tities  defined  herein  are  ex^re^sed  in  the  International  System  of  Units  (SI);  how- 


ever,  some  measurements  were  made  in  U.S.  Customary  Units  and  were  converted  by  fac 
tors  given  in  reference  7. 
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CT  thrust  coefficient,  Ttirust/puR  (21^) 

f  LV  burst  signal  frequency,  Hz 

f0R  Bragg  shift  frequency,  Hz 

Lf  fringe  spacing,  m 

number  of  velocity  measurements  in  ith  histogram  interval  as  percent  of 
total  number  of  measurements 
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order  number;  harmonic  of  one  rotor  revolution 
rotor  radius  (0.5  m) 
rotor-fuselage  spacing,  m 

local  velocity  sample  vector  aligned  with  x-axis,  m/sec  (see  fig.  1) 
ensemble  mean  of  u  measurement  component,  m/sec 
velocity  fluctuation  in  u  measurement  component,  m/sec 
local  velocity  sample  vector  aligned  with  z-axis,  m/sec  (see  fig.  1) 
ensemble  mean  of  v  measurement  component,  m/sec 
velocity  fluctuation  in  v  measurement  component,  m/sec 
free-stream  velocity,  m/sec 
fuselage  width,  m 

coordinate  system  relative  to  rotor  hub,  m 
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autocovariance  function,  m  /sec 

time  between  arrivals  of  sample  measurements,  described  in  azimuth 
intervals,  deg 


cross-beam  angle,  deg 

beam  wavelength,  m 

helicopter  advance  ratio,  V  /2uQR 

00 

free-stream  density,  kg/m^ 
rotor  rotational  speed,  Hz 


.—yrruai  I— 1 "  “ 

ret* 

HlXS 

UtlC 

justtflcatlon- 


pistriiwt  ion/ 


Superscript: 


peak  value 


TEST  APPARATUS 
Laser  Velocimeter 

The  optical  system  used  for  this  investigation  was  a  dual-color  four-beam 
fringe-type  laser  velocimeter  (LV)  operating  in  the  backscatter  mode.  A  schematic  of 
this  system  is  shown  in  figure  2.  Viewed  along  the  optical  axis,  the  diamond-shaped 
beam  pattern  was  composed  of  two  vertically  spaced  blue  beams  at  488.0  x  1O-9  m 
wavelength  and  two  horizontally  spaced  green  beams  at  514.5  x  10-9  m  wavelength. 

These  two  wavelengths  were  selected  by  a  color  separator  from  the  multiwavelength 
output  of  an  argon  ion  laser  operating  at  a  total  output  power  of  3  w.  The  two  beams 
were  routed  into  a  standard  two-component  LV  transmission  package.  There,  each  beam 
was  split  into  two  beams  of  equal  intensity,  and  one  beam  from  each  pair  was  passed 
through  an  acousto-optic  modulator  (Bragg  cell)  .  The  beam  spacings  were  then  stepped 
down  from  50  mm  to  22  mm,  which  was  the  distance  required  for  the  four-hole  mirror 
and  initial  lens  of  the  unique  afocal  zoom-lens  system.  This  zoom-lens  system  was 
designed  to  maintain  a  constant  cross-beam  angle  at  the  sample  volume  throughout  the 
zoom  range  (ref.  8).  The  system  has  three  lenses:  two  positive  lenses  separated  by 
a  distance  that  is  the  sum  of  their  focal  lengths,  and  a  negative  scan  lens.  The 
final  objective  lens  is  250  mm  in  diameter,  and  the  negative  scan  lens  is  32  mm  in 
diameter. 

A  summary  of  the  optical  characteristics  of  the  LV  system  is  presented  in 
table  I.  The  lens  system  maintained  a  constant  collecting  solid  angle  in  this 
coaxial  configuration  to  provide  a  constant  signal  amplitude.  The  final  folding 
mirror  was  the  prime  translation  device  of  the  sample  volume  for  this  investigation 
and  was  capable  of  two  degrees  of  freedom  about  a  horizontal  and  a  vertical  axis. 

At  the  crossover  point  if  the  four-beam  intersection  in  the  flow,  an  ellipsoidal 
sample  volume  was  formed  which  contained  a  distinctive  fringe  pattern,  as  described 
in  reference  9.  The  internal  fringe  spacing  L^r  was  a  function  of  wavelength  \ 
and  angle  9  between  the  crossing  beams  at  the  sample  volume: 


Lfr  2  sin(9/2) 


As  a  particle  passed  through  the  sample  volume,  light  was  scattered  with  an 
intensity  proportional  to  the  incident  intensity  in  the  fringe  pattern.  Part  of  the 
scattered  light  was  collected  by  the  transmitting  optics.  The  collected  scattered 
light  was  deflected  to  the  receiver  package  at  the  four-hole  mirror  between  the  nega¬ 
tive  lens  of  the  zoom  system  and  the  transmitting  optics  package.  In  the  receiver 
package,  the  two  combined  wavelengths  were  separated  and  directed  by  a  dichroic 
mirror  to  two  separate  photomultiplier  tubes  (PMT's),  where  each  optical  signal  was 
converted  to  an  electronic  signal,  and  a  classical  signal  burst  was  generated.  (See 
ref.  9.) 


The  Bragg  cell  provided  directional  information  for  the  velocity  measurement. 
In  this  particular  test,  the  effective  Bragg  shift  was  electronically  reduced  from 
40  MHz  to  1  MHz  for  improved  resolution,  since  the  velocities  expected  from  this 


rotor  were  below  8  m/sec.  The  velocity  calculated  from  the  LV  measurements?  is  a 
direct  function  of  the  fringe  spacing  described  earlier,  the  effective  Bragg  shift 
frequency  fBR,  and  the  frequency  f  within  the  LV  burst,  as  measured  by  the  elec¬ 
tronic  counters: 


u  or  v  =  (f  -  fBR)Lfr  (2) 

This  LV  system  was  installed  just  outside  the  tunnel  flow  in  the  test  chamber  of 
the  Langley  4-  by  7-Meter  Tunnel,  which  was  configured  as  an  open-throat  tunnel  with 
the  floor  in  place.  Hie  beam-crossing  optics  had  a  focal  range  from  2.45  to  4.9  m 
and  a  collecting  solid  angle  of  0.00132  sr.  The  sample  volume  was  0.63  cm  long  and 
0.18  mm  in  diameter,  and  was  positioned  at  desired  locations  within  the  rotor  plane 
using  the  pan-tilt  folding  mirror  and  the  zoom-lens  system.  The  photograph  of  the  LV 
optics  platform  in  figure  3  shows  the  laser  beam  paths  and  the  rotor  model  in  the 
background. 

The  LV  electronics  system  was  essentially  identical  to  that  described  in  refer¬ 
ence  '.0,  with  the  exception  of  the  LV  buffer  interface.  The  buffer  interface  used  in 
these  tests  was  capable  of  storing  two  channels  of  velocity  data;  however,  instead  of 
storing  the  elapsed  time  between  data  arrivals,  the  new  buffer  captured  rotor  azimuth 
simultaneously  with  the  two-component  velocity  measurements.  The  buffer  interface 
required  that  the  two  velocity  components  oe  coincident  in  time  in  order  for  the 
buffer  to  accept  the  measurement  as  valid  and  to  sample  the  azimuth  position.  Opera¬ 
tion  of  the  LV  buffer  interface  was  similar  to  that  described  in  reference  10,  but  in 
this  experiment  30  separate  buffer  cycles  were  initiated  at  each  desired  measurement 
position  instead  of  just  1.  Each  computer-controlled  cycle  was  halted  either  when 
the  interface  memory  was  full  (4096  data  points)  or  after  1  minute  of  elapsed  time. 

Particles  were  provided  by  the  smoke  generator  system  described  in  reference  10. 
The  particle  size  distribution  is  given  in  figure  4  (fig.  8(a)  of  ref.  10).  It  is 
important  to  note,  however,  that  the  optics  system  for  this  experiment  was  totally 
different  from  that  described  in  reference  10,  and  therefore  the  overall  measurement 
probability  distributions  s»h''**o  in  figures  8(b)  and  8(c)  of  reference  10  do  not  apply 
in  this  case. 


Rotor  Model 

The  model  used  in  this  investigation  was  an  externally  powered  1-m-diameter 
four-blade  rotor  with  three  fuselages  of  various  widths.  The  rotor  was  driven  by  an 
external  electric  motor  mounted  on  the  model  support  system  by  a  *?ix-component  bal¬ 
ance.  (See  fig.  1(b).)  A  fairing  was  provided  over  the  motor  and  balance  to  mini¬ 

mize  flow  interference  caused  by  mounting  hardware.  (See  fig.  3.)  Three  fuselages 
were  provided  (fig.  5),  all  1.0  m  long,  with  maximum  fuselage  widths  of  6.25  cm 
(0.125R),  12.5  cm  (0.250R),  and  18.75  cm  (0.375R)  (fig.  6).  In  this  configuration, 
the  rotor-fuselage  spacing  and  fuselage  width  could  be  varied  to  investigate  the 
effects  of  these  parameters  on  the  rotor  downwash  flow  field.  The  separation  dis¬ 
tances  and  fuselage  widths  used  for  this  investigation  are  shown  in  figure  6. 

The  rotor  system  was  operated  at  a  thrust  coefficient  CT  of  0.006,  a  rotor 
rotational  speed  Q  of  25  revolutions  per  second,  and  an  advance  ratio  ^  of  0.053, 

with  the  rotor  tip  path  plane  parallel  to  the  tunnel  horizontal  reference.  LV  meas¬ 

urements  were  obtained  in  a  plane  that  was  parallel  to  the  tip  path  plane  and  inter- 


sected  the  hub  and  flapping  hinges.  Since  the  model  was  mounted  upside  down,  the 
rotor  blades  coned  downward  due  to  thrust.  Therefore,  measurements  between  rotor 
blades  had  to  be  obtained  by  moving  the  LV  sample  volume  downward  from  the  rotor  hub 
plane.  Rotor  azimuth  was  provided  by  enhancing  a  72-per-revolution  (1  +  71)  pulse 
train  to  obtain  digital  levels  varying  from  0  to  720.  These  levels  were  then  asso¬ 
ciated  with  0,5°  rotor  azimuth  increments  from  0°  to  360°.  Rotor  performance  was 
monitored  continually  and  was  maintained  as  closely  as  possible  to  test  conditions 
throughout  each  data  point.  Selected  rotor  and  tunnel  parameters  were  passed  from 
the  existing  tunnel  data  acquisition  computer  to  the  LV  data  acquisition  computer 
prior  to  each  LV  buffer  cycle. 


DATA  ACQUISITION  AND  REDUCTION 
Statistical  Quantities 

The  LV  samples  a  velocity  field  in  time  with  the  variability  of  a  Poisson  dis¬ 
tribution.  (See  ref.  1.)  The  velocity  measurement  from  a  particle  is  only  one  sam¬ 
ple  of  velocity  within  the  probability  density  function  of  the  flow  field.  If  the 
standard  deviation  of  the  velocity  measured  is  large,  the  probability  of  this  single 
measurement  being  a  good  representation  of  the  expected  value  of  the  true  flow  veloc¬ 
ity  is  statistically  rather  low.  Therefore,  it  is  necessary  to  acquire  an  adequate 
number  of  samples  in  order  to  completely  define  the  statistical  probability  density 
function  of  the  velocity  measurements.  During  the  data  reduction  process,  two 
assumptions  are  made:  (1)  the  particles  embedded  in  the  flow  are  not  only  randomly 
dispersed  in  space  but  are  also  randomly  dispersed  in  the  velocity  field;  and  (2)  the 
measurement  sample  ensemble  obtained  over  a  finite  period  of  time  is  a  good  statis¬ 
tical  representation  of  the  stationary-flow  condition  at  the  measurement  location. 

In  this  investigation,  it  will  be  shown  that  flow  phenomena  affecting  this  stationar- 
ity  can  have  serious  effects  on  the  velocity  measurements  obtained. 

This  investigation  was  conducted  to  determine  the  periodic  variation  and  stan¬ 
dard  deviation  of  local  velocity  with  rotor  azimuth.  These  azimuth-dependent  quan¬ 
tities  could  be  obtained  by  sorting  the  velocity  information  into  ensembles  of  equal 
azimuth  (0.5°  resolution)  and  performing  the  statistical  calculations  described  in 
reference  10  on  each  azimuth  ensemble.  Time-averaged  mean  velocity  calculations  were 
made  based  on  the  simple  arithmetic  mean, 


Zv, 

- _ i _ 

Number  of  velocity  measurements 

It  was  not  necessary  to  make  corrections  for  velocity  bias  or  Bragg  cell  bias,  as 
noted  in  reference  10. 


Error  Analysis 

The  overall  measurement  precision  was  obtained  by  determining  the  accuracy  of 
each  variable  in  the  system  that  could  affect  the  accuracy  of  the  velocity  measure¬ 
ments.  Reference  11  provides  a  complete  description  of  the  sources  of  errors 
involved  in  laser  velocimetry,  as  well  as  of  the  error  analysis  method.  These  error 
sources,  presented  in  table  II,  yield  an  overall  effective  total  bias  error  in  the 
mean  from  -1 .0  to  +0.74  percent.  These  values  were  calculated  from  an  algebraic  sum 


of  the  partial  bias  errors.  The  total  effect  of  random  error  on  the  mean  was 
±0.58  percent  uncertainty.  This  value  was  obtained  by  taking  the  square  root  of  the 
sum  of  the  squates  of  the  partial  random  errors. 

Instrument-related  measurement  uncertainties  were  caused  by  cross -beam-angle 
uncertainty,  nonparallelism  of  fringes,  clock  synchronization,  and  quantizing  error. 
The  cross-beam-angle  uncertainty,  which  resulted  from  inaccuracies  in  the  measurement 
of  the  cross-beam  angle,  produced  a  constant  but  unknown  bias  error  in  the  measure¬ 
ment  of  particle  velocity.  However,  these  inaccuracies  did  not  produce  any  random 
uncertainties  in  the  mean  velocity  measurements.  Hie  uncertainty  due  to  nonparallel¬ 
ism  of  fringes  yielded  a  bias  and  a  random  uncertainty  in  the  mean  velocity  measure¬ 
ment.  Hie  synchronization  between  the  signal  burst  and  the  clock  in  the  high-speed 
burst  counter  caused  a  bias  error  and  a  negligible  random  error  in  the  mean  velocity 
measurement.  For  this  test,  the  effects  of  the  quantizing  error  caused  by  the  digi¬ 
tization  of  time  in  the  high-speed  burst  counter  were  negligible  for  the  measurement 
of  mean  velocity. 

In  a  large  velocity  gradient,  velocity  measurement  errors  can  occur  due  to  the 
length  of  the  sample  volume,  which  allows  contributions  from  the  high-velocity  flow 
at  one  end  and  the  low-velocity  flow  at  the  other  end  to  be  included  in  the  measure¬ 
ment.  If  the  gradient  is  linear  over  the  length  of  the  sample  volume,  as  is  assumed 
in  the  present  study,  the  mean  velocity  measurement  is  not  affected.  Further,  the 
results  of  this  test  show  that  the  flow  was  below  10  m/sec  with  a  turbulence  intens¬ 
ity  of  less  than  10  percent,  indicating  that  the  velocity  gradients  were  very  small. 
Thus,  in  general,  the  error  caused  by  velocity  gradients  can  be  neglected. 

Since  the  LV  measures  particle  velocities  and  not  the  gas  velocity,  the  final 
measurement  accuracy  is  dependent  on  the  ability  of  the  particle  to  follow  the  flow 
faithfully.  In  the  present  study,  kerosene  droplets  were  used  as  the  seed  material. 
The  particle  size  distribution  is  shown  in  figure  4.  Based  on  the  results  reported 
in  reference  11,  these  particles  will  follow  a  velocity  gradient  of  1500  (m/sec)/sec 
within  approximately  5  percent.  Thus,  it  is  concluded  that  these  particles  will 
faithfully  follow  the  flow  field  studied  in  this  test. 


TEST  PROCEDURES 

This  investigation  was  conducted  in  the  Langley  4-  by  7-Meter  Tunnel.  Rotor 
wake  velocity  measurements  were  obtained  in  at  least  8  locations  in  the  plane  of  the 
hub  for  10  different  rotor-fuselage  configurations:  rotor  only,  and  3  rotor-fuselage 
spacings  (S/R  =  0.224,  0.275,  and  0.324)  each  for  3  different  fuselage  widths 
(W/R  =-  0.125,  0.250,  and  0.375).  (See  table  AI  for  a  listing  of  measurement  loca¬ 
tions.)  Most  of  these  measurements  were  obtained  along  the  x-axis  upstream  of  the 
rotor  hub.  Some  data  were  acquired  very  near  and  in  the  surface  of  the  cone  inscribed 
by  the  rotor  blades;  this  required  making  velocity  measurements  between  rotor  blades. 

The  LV  data  were  acquired  as  particles  passed  through  the  sample  volume.  Par¬ 
ticles  were  injected  into  the  flow  at  sufficient  density  and  size  to  produce  a  high 
arrival  rate.  As  the  counter-processor  validated  each  measurement,  it  was  accepted 
and  stored  in  the  high-speed  buffer  along  with  the  appropriate  rotor  azimuth.  A 
coincidence  time  restriction  of  1.0  x  10~®  sec  was  placed  on  the  u-  and  v-component 
measurements.  This  computer-controlled  process  continued  until  the  buffer  was  filled 
or  the  time  limit  was  exceeded.  The  entire  measurement  sequence  was  repeated 
30  times  for  each  measurement  location. 


RESULTS  AND  DISCUSSION 


Examples  of  velocity  measurements  obtained  for  each  sample  volume  location  are 
presented  in  the  appendix.  Only  selected  samples  of  these  data  will  be  discussed  in 
the  text,  along  with  simulated  samples  of  predicted  time-dependent  velocity  profiles. 
Power  spectral  density  and  autocorrelation  function  estimates  of  the  experimental  and 
simulated  data  sets  are  presented  and  discussed.  Time-averaged  calculations  of  the 
experimental  data  acquired  are  presented  in  figures  7  and  8. 


Effects  of  Rotor-Fuselage  Configuration  on  Time-Averaged 
Velocity  Measurements 

Data  presented  in  figures  7  and  8  are  time-averaged  samples  of  data  obtained 
along  the  x-axis  in  the  plane  of  the  rotor  hub  (i.e.,  the  measurements  were  made 
between  the  rotor  blades  and  the  fuselage).  (See  fig.  1.)  These  data  are  presented 
as  horizontal  (u)  and  vertical  (v)  velocity  components  as  a  function  of  distance  from 
the  hub  normalized  by  rotor  radius  R.  The  effect  of  body  width  at  a  fixed  rotor- 
fuselage  spacing  is  presented  in  figure  7.  The  primary  effect  of  body  width  occurs 
at  the  inner  radius  of  the  rotor,  which  is  the  area  most  affected  by  the  fuselage 
presence.  As  expected,  the  vertical  velocity  component  decreased  and  the  horizontal 
velocity  component  increased  with  increasing  fuselage  width.  This  had  a  channeling 
effect  on  the  flow. 

The  effects  of  rotor-fuselage  separation  for  each  fuselage  width  are  presented 
in  figure  8.  The  presence  of  the  0.1 25R  body  at  different  distances  from  the  rotor 
is  barely  evident  in  the  velocity  measurements  obtained  (fig.  8(a)),  but  the  effects 
of  separation  distance  for  the  other  two  bodies  are  larger  (figs.  8(b)  and  (c)).  In 
general,  the  rotor  wake  is  being  channeled  between  the  fuselage  and  the  rotor  in  the 
close-coupled  configuration  (S/R  =  0.224),  as  evidenced  by  the  higher  horizontal  and 
lower  vertical  velocity  components  at  x/R  =  0,1.  (See  fig.  8(b).)  Slightly 
increased  vertical  and  decreased  horizontal  velocity  components  for  the  close-coupled 
rotor-fuselage  configuration  at  x/R  =0.3  indicate  acceleration  of  the  rotor  wake 
conditions  around  the  nose  of  the  fuselage.  These  effects  can  be  observed  in  the 
data  for  the  0.250R  and  0.375R  bodies;  the  wider  body  obviously  shows  larger  magni¬ 
tude  effects. 


Analysis  of  Time-Dependent  Data 

It  appears  from  the  LV  data  presented  in  the  appendix  that  no  real  periodic 
velocity  variation  with  rotor  blade  passage  (azimuth)  can  be  defined.  Since  the 
prime  objective  of  this  program  was  to  identify  this  time-dependent  variation,  it 
would  seem  that  the  project  was  not  successful.  However,  to  fully  understand  the  LV 
process  and  its  application  to  rotor  wake  measurements,  it  is  imperative  to  identify 
the  cause  of  the  apparent  lack  of  time -dependent  variation  in  the  data.  Three  possi¬ 
bilities  can  explain  this  result:  (1)  the  azimuth  indicator  was  not  operating  cor¬ 
rectly  and  provided  erroneous  azimuth  measurements;  (2)  the  velocity  truly  did  not 
vary  with  time;  or  (3)  some  flow  phenomenon  caused  variability  in  the  velocity  time 
history  during  the  measurement  period,  thus  washing  out  the  time-dependent  velocity 
measurements  (i.e.,  there  was  a  lack  of  stationarity  in  the  flow  field). 

To  address  the  first  possibility,  a  sample  velocity  measurement  was  selected  for 
analysis.  Figure  9  shows  a  sample  (one  buffer  load)  of  the  velocity  measurements 
obtained  when  the  sample  volume  was  located  in  the  path  of  the  rotor  blades.  Veloc- 


ity  components  as  a  function  of  azimuth  are  presented  in  figure  9(a),  and  figure  9(b) 
presents  the  velocity  histogram  (or  probability  density  function)  of  these  compo¬ 
nents.  Each  dot  in  figure  9(a)  represents  a  single  velocity  measurement.  Note  the 
apparent  voids  in  the  velocity  data  at  azimuth  angles  from  30°  to  60°,  120°  to  140°, 
210°  to  230°,  and  300°  to  320°.  The  voids  are  caused  by  the  passage  of  each  blade 
through  the  sample  volume.  If  the  azimuth  measurement  had  been  malfunctioning,  these 
voids  would  not  exist  here  or  in  the  other  data  acquired  during  the  time  the  beams 
were  being  clipped  by  the  rotor  blades.  Note  also  that  although  velocity  measure¬ 
ments  were  made  between  the  blades,  there  are  still  no  apparent  time -dependent  veloc¬ 
ity  characteristics.  It  is  evident  that  the  azimuth  measurement  was  functioning 
properly  and  is  not  a  probable  cause  of  the  lack  of  time-dependent  velocity  varia¬ 
tions  at  the  selected  measurement  locations. 

The  second  possibility,  that  the  velocity  did  not  vary  with  time,  was  examined 
by  calculating  the  wake  velocities  using  an  analytical  prediction  method  frequently 
employed  in  rotor-wake  analyses.  This  method,  developed  by  Landgrebe  (ref.  12),  was 
used  to  estimate  the  time-dependent  velocity  characteristics  for  this  rotor  wake 
without  fuselage.  The  classical  wake  calculations  (undistorted  helix)  used  for  this 
particular  program  did  not  allow  distortion  of  the  wake  shape  due  to  wake  interac¬ 
tion.  Precise  correlation  of  the  velocity  magnitude  with  the  LV  measurement  was  not 
attempted  in  this  case,  since  only  an  example  of  the  time-dependent  nature  of  the 
rotor  wake  at  the  desired  location  was  required.  Thus,  the  velocity  magnitudes  did 
not  have  to  be  very  accurate  for  this  prediction. 

Three  velocity  predictions  are  presented  in  figure  10.  The  solid  line  in  fig¬ 
ure  10(a)  represents  the  estimated  velocity  time  history.  For  this  example,  the 
simulation  calculation  assumed  a  normal  density  function  distribution  for  the  veloc¬ 
ity  ensemble  at  each  azimuth,  with  a  standard  deviation  of  1  percent  in  both  the 
vertical  and  the  horizontal  component.  The  dots  in  this  figure  were  simulated  using 
a  Poisson  time  distribution  (ref.  1)  and  represent  samples  of  the  assumed  velocity 
distribution.  The  waveform  character  of  the  velocity  time  histories  is  preserved, 
probably  because  the  turbulence  levels  were  random  fluctuations.  Random  turbulence, 
however,  can  be  a  serious  matter  if  the  intensities  are  large,  as  seen  in  the  simu¬ 
lated  calculations  obtained  with  an  assumed  turbulence  of  10  percent  (fig.  10(b)). 

In  this  case,  no  velocity  waveform  characteristic  would  be  apparent  unless  it  were 
already  known.  It  should  be  noted  that  actual  turbulence  intensity  level  in  the 
closed  test  section  of  the  4-  by  7-m  tunnel  is  0.3  percent  for  the  horizontal  com¬ 
ponent  and  1 .0  percent  for  the  vertical  component,  which  is  well  within  the  limits 
for  retaining  the  velocity  waveform  characteristics. 

For  this  particular  investigation,  the  tunnel  was  operated  with  the  test  section 
open.  After  the  investigation,  measurements  of  tunnel  turbulence  values  indicated 
that  a  pulsation,  or  discrete-frequency  velocity  fluctuation,  was  present  at  about 
1.5  Hz.  The  magnitude  of  this  discrete-frequency  pulsation  at  the  lowest  speed 
tested  was  on  the  order  of  8  percent  in  the  vertical  component  and  5  percent  in  the 
horizontal  component  with  reference  to  the  free-stream  velocity.  (See  fig.  11.) 
Unfortunately,  turbulence  measurements  were  not  obtained  at  velocities  as  low  as 
those  used  in  this  investigation. 

To  consider  a  worst-case  possibility,  the  10-percent  turbulence  intensity  was 
chosen  for  the  open-tunnel  simulation  calculations.  Such  large  vertical  and  horizon¬ 
tal  discrete-frequency  fluctuations  could  result  in  a  large  variation  (as  much  as  4°) 
in  the  angle  of  the  flow  presented  to  the  rotor  system,  assuming  that  the  velocity 
components  were  shifted  in  phase  180°  relative  to  the  pulsation.  The  Landgrebe  pre¬ 
diction  technique  (ref.  12)  was  used  to  estimate  the  periodic  variation  in  the 
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azimuth-dependent  velocity  at  one  point  in  the  rotor  wake.  Based  on  these  predic¬ 
tions,  simulated  LV  measurements  were  generated  for  a  model  in  the  open  test  section. 
These  simulated  measurements  were  obtained  using  the  same  calculations  described  pre¬ 
viously,  but  with  the  addition  of  periodic  {1.5  Hz)  free-stream  velocity  fluctua¬ 
tions.  As  can  be  seen  from  figure  10(c),  no  time-dependent  velocity  characteristic 
is  evident. 

These  results  indicate  that  the  periodic  free-stream  velocity  fluctuations 
caused  a  nonstationary  flow  condition  at  the  measurement  location.  The  LV  (or,  for 
that  matter,  any  instrument  physically  located  at  a  point  in  this  environment)  would 
have  measured  the  velocity  characteristics  of  this  nonstationary  oscillating  rotor 
wake  in  addition  to  those  of  the  stationary  wake  desired.  Measurements  obtained 
experimentally  at  the  same  locations  used  in  calculating  the  simulated  results  shown 
in  figure  10  (i.e.,  x/R  =  0.3,  y/R  =  0,  and  z/R  =  0)  also  show  this  lack  of  time- 
dependent  velocity  characteristic.  (See  fig.  12.)  The  remarkable  similarity  of  the 
data  scatter  in  these  figures  leads  to  the  postulation  that  the  high  turbulence  level 
and  its  periodic  nature  obscured  the  time-dependent  characteristics  of  the  velocity 
measurements.  However,  it  is  impoitant  to  note  that  this  is  only  a  definite 
possibility  and  has  not  yet  been  proven. 


Spectral  Characteristics 

Tfc>  further  analyze  the  effects  of  this  periodic  turbulence,  the  data  (both 
experimental  and  simulated)  were  transformed  into  the  frequency  domain  using  autoco¬ 
variance  and  fast  Fourier  transform  ( FFT)  techniques.  It  is  well  known  that  an  auto¬ 
covariance  of  a  random-noise  data  set  that  has  a  periodic  component  buried  within  it 
will  provide  a  measure  (frequency  and  amplitude)  of  the  periodic  component  expressed 
by  the  autocovariance  function.  An  FFT  of  this  autocovariance  function  will  then 
give  an  estimate  of  the  power  spectral  density  of  the  data  set.  It  has  been  demon¬ 
strated  (ref.  13)  that  direct  calculation  of  the  autocovariance  function  is  necessary 
for  LV  data  transformation. 

The  autocovariance  function  reveals  any  coherent  relationship  of  the  data  under 
investigation  as  a  function  of  time.  In  this  test,  the  autocovariance  function  was 
determined  as  a  function  of  azimuth,  so  that  minor  discrepancies  due  to  small  varia¬ 
tions  in  rotor  speed  would  not  affect  the  calculations.  In  laser  velocimetry,  the 
individual  measurements  of  velocity  occur  at  random  times;  this  requires  a  knowledge 
of  the  time  (in  this  case,  azimuth  angle  A4>)  between  successive  measurements.  All 
possible  velocity  pairings  were  made  for  a  total  differential  azimuth  angle  of  less 
than  one  revolution.  The  cross  products  were  then  summed,  as  was  the  number  of  times 
each  differential  azimuth  angle  occurred.  Normalization  of  the  cross-product  summa¬ 
tions  by  the  number  of  occurrences  of  each  angle  yielded  the  final  autocovariance 
function.  This  process  may  be  expressed  as 
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where 


N  number  of  measurements  (i.e.,  particles) 

u(j)  jth  velocity  measurement  (with  mean  removed) 

u(j+k)  velocity  measurement  kA(Jj  from  j  (with  mean  removed) 

I(j)  binary  indicator  (I(j)  =  1) 

I(j+k)  binary  indicator  (I(j+k)  =1  if  there  is  a  velocity  measurement  at  j  +  k, 

and  I(j+k)  =  0  otherwise) 

A<|j  autocovariance  delay  angle 

Once  the  autocovariance  function  was  calculated,  it  was  multiplied  by  a  Bartlett 
window  function  to  reduce  the  effect  of  the  side  lobes  in  the  spectral  window.  FFT 
techniques  were  applied  to  the  results  to  yield  the  power  spectral  density  (PSD)  of 
the  flow  field  under  observation. 

The  simulated  data  set  presented  in  figure  10(a),  which  represented  a  sample 
condition  in  the  closed  tunnel,  was  transformed  in  this  manner.  The  results  are 
presented  in  figure  13.  These  simulated  data  can  be  seen  to  have  a  very  strong  peri¬ 
odic  content  in  the  time  domain  (fig.  10(a))  -  4  per  revolution  in  the  u  component 
and  8  per  revolution  in  the  v  component.  The  autocovariance  function  in  fig¬ 
ure  13(a)  shows  the  classical  sinusoidal  variation  of  the  discrete  frequency  in  the 
u  and  v  components.  The  estimates  of  power  spectral  density  that  were  produced  by 
fast  Fourier  transforms  of  the  autocovariance  function  are  presented  in  figure  13(b). 
These  results  are  presented  as  a  function  of  order  number  (or  multiples  of  rotor 
revolutions).  The  PSD  shows  the  discrete-frequency  velocity  characteristics  at  those 
values  observed  in  the  velocity  time  histories.  (See  fig.  10(a).) 

By  performing  this  same  transformation  on  the  data  in  figure  10(c),  any  dis¬ 
crete-frequency  velocity  components  that  exist  will  be  seen  ir.  the  autocovariance 
function  and  the  PSD.  This  data  set  simulated  data  acquisition  in  the  tunnel  with 
the  test  section  open  at  a  point  where  periodic  turbulence  washed  out  any  apparent 
time-dependent  velocity  characteristics.  The  results  of  spectral  transformation  are 
presented  in  figure  14.  No  apparent  periodic  content  is  evident  in  the  autocovari¬ 
ance  function,  but  a  very  weak  8-per-revolution  frequency  content  may  be  present,  as 
shown  in  the  v  component  in  figure  14(b). 

The  results  of  these  calculations  on  the  simulated  LV  data  suggest  two  important 
conclusions.  First,  if  a  discrete-frequency  content  exists  in  the  random  LV  data, 
the  spectral  transformation  process  described  should  identif}  it.  Second,  if  the 
apparent  lack  of  time-dependent  velocity  characteristics  in  the  experimental  data  is 
in  fact  a  result  of  the  periodic  turbulence  in  the  open  test  section,  then  spectral 
transformation  of  the  experimental  data  will  show  no  periodic  content.  This  trans¬ 
formation  was  performed  on  the  experimental  LV  data,  and  the  results  are  presented  in 
figure  15.  No  periodic  content  is  observable  in  either  the  autocovariance  function 
or  the  power  spectral  density.  Data  for  many  other  rotor  configurations  and  LV  mea¬ 
surement  positions  were  transformed  into  the  frequency  domain  with  similar  results. 


CONCLUDING  REMARKS 


An  experimental  investigation  was  conducted  using  a  laser  velocimeter  (LV)  to 
determine  the  effects  of  various  fuselage  widths  and  rotor-fuselaqe  separations  on 
time-averaged  and  detailed  time-dependent  rotor  wake  velocity  characteristics.  The 
results  of  the  investigation  were  ambiguous  in  that  no  apparent  time-dependent  veloc¬ 
ity  characteristics  were  observed  in  the  experimental  data.  Simulated  LV  sampling  of 
a  theoretically  determined  velocity  time  history  demonstrated  the  probable  cause  of 
these  results.  Using  a  recognized  analytical  prediction  method,  it  appeared  that 
velocity  variations  should  have  existed  at  the  measurement-  locations  chosen.  Time- 
dependent  variations  in  velocity  should  have  been  preserved  for  turbulence  charac¬ 
teristics  typical  of  the  closed  test  section.  However,  these  time-dependent  velocity 
characteristics  were  totally  washed  out  in  the  open-throat  tunnel  configuration. 
Autocovariance  calculations  and  power  spectral  density  estimates  were  used  to  verify 
that  the  time-dependent  velocity  characteristics  were  not  merely  buried  within  the 
observed  random- velocity  time  history. 

Based  on  these  results,  it  is  most  probable  that  the  discrete-f requency  pulsa¬ 
tion  existing  in  the  facility  in  the  open-throat  configuration  at  the  speeds  tested 
set  up  a  nonstationary  flow  field,  thus  presenting  different  test  conditions  to  the 
LV  instrument  over  the  time  period  of  an  ensemble  measurement.  When  this  program  was 
planned,  flow  stationarity  in  the  test  facility  was  assumed,  and  flow  fluctuations  at 
frequencies  below  the  rotor  rotational  frequency  were  ignored.  The  data  acquisition 
technique  was  not  designed  to  record  elapsed  time  between  measurements  larger  than 
the  time  for  one  revolution.  This  limitation  made  it  impossible  to  isolate  velocity 
fluctuations  at  frequencies  less  than  the  rotor  rotational  frequency.  Since  these 
velocity  fluctuations  could  not  be  separated  from  those  at  frequencies  of  interest, 
they  could  not  be  prevented  from  contributing  to  the  scatter  and  the  amplitude  uncer¬ 
tainty  in  the  rotor-related  periodic  content.  However,  it  is  important  to  point  out 
that  operation  at  higher  rotor  disk  loadings  or  different  tunnel  velocities,  where 
these  periodic  turbulence  levels  are  lower,  might  have  produced  better  results.  It 
is  recommended  that  any  future  programs  for  which  large  time-dependent  velocity  gra¬ 
dients  exist  as  a  stationary  process  should  be  conducted  in  the  closed  test  section 
of  the  facility.  Programs  using  the  open  test  section  should  be  examined  closely  to 
insure  that  this  discrete-frequency  turbulence  will  not  create  a  nonstationary  flow 
field  that  may  render  the  measurements  useless. 
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TABLE  I.-  OPTICAL  CHARACTERISTICS  OF  LV  SYSTEM 


Input  beam  diameter,  mm  . . . . . . .  1.8 

Input  beam  separation,  mm  . . .  22 

Output  beam  diameter,  mm  . . .  8.1  to  17.3 

Output  beam  separation,  cm . . .  9.7  to  20.3 

Included  angle,  0,  deg  . . .  2.36 

Scan  range,  m  . . .  2.45  to  4.9 

Diameter  of  objective  lens,  mm  . . . . .  250 

Approximate  number  of  fringes  in  sample  volume  .  15 

Fringe  spacing,  m: 

X  =  488.0  x  10~9  m  .  11.85  x  10-6 

X  =  514.5  x  10-9  m  .  12.5  x  10-6 

Sample  volume: 

Length,  mm  . . .  6.3 

Diameter,  mm  . . .  0.18 


Collection  solid  angle,  sr 


0.00132 


TABLE  II.-  MEAN  VELOCITY  ERROR  SOURCES 


Error  source 

Bias,  percent 

Random, 

percent 

Velocity  bias 

Corrected 

NA 

Bragg  bias 

Corrected 

NA 

Cross-beam-angle  uncertainty 

±0.87 

NA 

Nonparallelism  of  fringes 

-0.28 

±0.42 

Time  jitter 

0 

0 

Clock  synchronization 

.15 

NA 

Quantizing  error 

(a) 

±.0001 

Velocity  gradient 

(a) 

(a) 

Measurement  location  uncertainty 

(b) 

(b) 

Particle  lag 

(a) 

(a) 

Statistical  uncertainty  (typical) 

NA 

±0.40 

Total  error 

+0.74  to  -1 .0 

±0.58 

(a)  Top  view. 


(b)  Side  view. 


Schematic  of  rotor  system  tested,  including  axis  system 


Schematic  of  laser  velocimeter  optics  package 


7-Meter  Tunnel. 
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Figure  7.-  Effect  of  body  width  on  velocities  measured 

S/R  =  0.224. 


Effect  of  rotor-body  spacing 


( rotor  only) 


Figure  8.-  Continued 


(b)  Closed  tunnel,  10  percent  random  tunnel  turbulence 


Figure  11.-  Turbulence  intensity  characteristics  of  open  test  section  of 

Langley  4-  by  7-Meter  Tunnel. 
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Velocity,  m/sec 


(a)  Velocity  components  as  a  (b)  Histograms  of  velocity  components 

function  of  azimuth. 

Figure  12.-  Laser  velocimeter  measurements  in  rotor-only  wake.  x/R  =  0.3; 
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(a)  Autocovariance  function. 


(b)  Power  spectral  density. 

Figure  13.-  Spectral  transformation  of  simulated  laser  velocimeter  values 
obtained  using  theoretically  determined  velocity  time  history.  Closed 
tunnel;  random  tunnel  turbulence. 


(a)  Autocovariance  function 


Figure  1 


(b)  Power  spectral  density. 

5.-  Spectral  transformation  of  laser  velocimeter  measurements  in 
rotor-only  wake.  x/R  =  0.3;  y/R  =  0;  z/R  =  0. 


APPENDIX 


SAMPLE  VELOCITY  MEASUREMENTS 

Samples  of  each  velocity  measurement  obtained  in  this  investigation  are  given  in 
this  appendix.  Each  figure  represents  one  buffer  load  (out  of  30)  for  each  rotor 
configuration  and  position  tested.  Each  dot  in  the  plots  represents  a  single  mea¬ 
surement,  and  since  the  data  system  was  restricted  to  coincidence  (data  accepted  only 
when  obtained  by  both  u  and  v  counters  simultaneously),  each  v  component  dot 
has  a  corresponding  u  component  dot.  Histograms,  or  probability  density  functions, 
are  presented  for  each  set  of  u-  and  v-component  measurements.  A  sample  interpreta¬ 
tion  of  these  histograms  can  be  found  in  appendix  B  of  reference  10. 

The  presentation  of  the  figures  for  the  data  acguired  in  this  investigation  is 
given  in  table  AI. 


Figure  A1 LV  data  for  rotor-only  test 
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Figure  A1 Continued 
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Figure  A2.-  LV  data  for  rotor  with  0.1 25R  body 
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Figure  A3.-  LV  data  for  rotor  with  0.250R  body 
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Figure  A3.-  Continued 
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Figure  A4.-  Continued 
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